Entrainment parameters in cold superfluid neutron star core 



Nicolas Chamel 
Copernicus Astronomical Center, Polish Academy of Science, 
ul. Bartycka 18, 00-716 Warszawa, Poland and 
LUTH, Paris Observatory, 5 place Jules Janssen, 92195 Meudon, France 



o 
o 

(N 

(D 

m 



Pawel Haensel 

Copernicus Astronomical Center, Polish Academy of Science, ul. 

(Dated: February 7, 2008) 



Bartycka 18, 00-716 Warsaw, Poland 



Hydrodynamical simulations of neutron star cores, based on a two fluid description in terms of 
a neutron-proton superfluid mixture, require the knowledge of the Andreev-Bashkin entrainment 
matrix which relates the momentum of one constituent to the currents of both constituents. This 
matrix is derived for arbitrary nuclear asymmetry at zero temperature and in the limits of small 
relative currents in the framework of the energy density functional theory. The Skyrme energy 
density functional is considered as a particular case. General analytic formulae for the entrainment 
parameters and various corresponding effective masses are obtained. These formulae are applied to 
the liquid core of a neutron star, composed of an homogeneous plasma of nucleons, electrons and 
possibly muons in beta equilibrium. 
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I. INTRODUCTION 

In a standard model of neutron star core, matter is a 
uniform plasma consisting of neutrons, of number den- 
sity n„, and a small admixture of protons and electrons, 
of number densities rip and rie, respectively. Electrons 
ensure the overall stability of the star by the condition 
of electroneutrality, n,, — rip, but play a negligible role in 
the mass transport because their mass is very small com- 
pared to the nucleon mass. If the electron Fermi energy 
exceeds the muon rest mass, muons are present in matter, 
but their density is always smaller than that of electrons. 
The electroneutrality condition is then rie -\- = Up. 
The electrically charged particles are strongly coupled to 
the magnetically braked solid crust. This builds up a 
lag between the neutrons and the protons which is only 
restored through glitch events. Neutron star cores are 
therefore described within two fluid models in terms of 
neutron and proton components which are superfluid in 
some density range when the temperature of the star falls 
below the corresponding critical temperatures. 

As a result of the nucleon- nucleon interactions, the mo- 
mentum TTq of each nucleon is not simply given by the 
corresponding velocity Vq (we use the convention that 
q = n,p for neutron, proton respectively) times the mass 
m (in the following we shall neglect the small mass differ- 
ence between neutrons and protons) but in general is a 
linear combination of the neutron and proton velocities. 
This is the so called Andreev-Bashkin entrainment effect 
0. The discussions in the literature have been usually 
obscured by the confusion between momentum and veloc- 
ity. Traditionally one introduces "superfluid velocities" 
Vq defined by 



in which it is recalled that the momenta TTq are defined by 
the partial derivative with respect to the nucleon current 
nnVn of the Lagrangian density Ajn^, UqVq} of the system 
[2 y| ■ The mass current of some given nucleon species q 



Pq = PqVq , Pq 



is then expressible as 



(2) 



(3) 
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in which pqqi is the (symmetric) entrainment or mass 
density matrix. Only one of these matrix elements has 
to be specified since the other elements can be obtained 
from the identities due to Galilean invariance 



Pnn ^ P np Pn ; Ppp ~^ Ppn — Pp ■ 



(4) 



Vq = TTq/m, 



(1) 



This matrix is a necessary ingredient in dynamical sim- 
ulations of neutron star cores, such as for instance the 
study of oscillation modes. The (static) equation of state 
and the entrainment matrix are usually obtained using 
different microscopic models. In earlier calculations and 
even recently, the mass density matrix is postulated to 
have some density dependence whose parameters are de- 
termined from rough estimates. 

Comer et al. have built a self-consistent equation 
of state in the framework of a minimal relativistic a — uj 
mean field model, ignoring non linear couplings between 
the meson fields which are however essential in order to 
reproduce nuclear properties such as the incompressibil- 
ity of nuclear matter. They have obtained semi analytical 
formulae for the entrainment parameters in the limit of 
small fluid velocities (compared to that of light), which 
even within this simple mean field model take a rather 
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complicated form. It is not clear that analytical formu- 
lae could still be obtained with realistic relativistic mean 
field models, taking into account self meson couplings 
and including as well the p meson which is required for 
a correct treatment of the symmetry energy. 

Despite the fact that non relativistic mean field mod- 
els have been widely applied in the study of terrestrial 
nuclei and in neutron stars, there has been no attempt 
so far to apply these models for the calculation of the 
mass density matrix. The purpose of the present work 
is therefore to fill this gap and to further investigate the 
density dependence of the various entrainment parame- 
ters and effective masses that have been introduced in 
the literature. 



II. ENTRAINMENT IN A MIXTURE OF 
FERMI LIQUIDS 

At zero temperature the entrainment effects have been 
shown to be independent of the nucleon pairing, giving 
rise to superfluidity 0, 0| . Even at finite temperatures 
well below the critical temperatures for the onset of su- 
perfluidity, pairing as well as thermal effects are very 
small 0. We can therefore ignore pairing interactions 
and restrict ourself to the limit of zero temperature. 

Borumand et al. (7j have shown how to obtain the 
entrainment matrix of a neutron-proton mixture in the 
framework of the Landau Fermi liquid theory. In what 
follows, we will limit ourselves to spin-unpolarized nu- 
clear matter. Therefore, spin indices will not appear in 
our formulae, and all quantities are to be understood as 
spin-averages. Under our assumptions, the change in the 
total energy density of the system due to a small current 
is expressed as 



i2n) 

/ ^ / ^f'i'{k,k'}Sn<^i^k}Sn('i'Hk'}, 



(27r)3 J (27r)3 



(5) 



it can be shown that the mass current = pgVq of each 
nucleon species is linearly related to both the neutron 
and proton superfluid velocities 



Pq ~ Pqq'Vq' , 



(7) 



where the (symmetric) entrainment matrix pqqi is given 
by 



qq> 



/3) 



(8) 



The (Landau) effective mass to® and the dimensionless 
Landau parameters J-f^ are defined respectively by 



1 de 
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m® h^ki^' dfc 



(9) 



(10) 



in which J\fq is the density of quasiparticle states at the 
Fermi surface. 



(9) 
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h 7T 



(11) 



and the parameters /^"^ are obtained from the Legendre 
expansion of the spin averaged quasiparticle interaction 



r'^'{k,k'} = Y,fr'pdco8e} 



(12) 



where 9 is the angle between the wave vectors k and k' 
lying on the corresponding Fermi surface. 

Alternative formulae of the entrainment matrix Q 
have been used in the literature, based on the decom- 
position of the Landau effective masses in the form 



TO + 6m. 



© 
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(5to„ 



(13) 



in which e^'^'jlc} is the energy of a quasiparticle {q = n,p 
for neutron and proton respectively) of wave vector A;, 
{k,k'} is the (spin averaged) interaction between the 
quasiparticles. Moreover, Sn^"^ ^{k'} denotes the change 
in the distribution function of quasiparticle states from 
that of the static (zero current) ground state charac- 
terized by the Heaviside functions Q{kp^ — k}, where 
kp^ is the Fermi momentum (in the units of h) kp^ = 
(37r^nq)^/3 jji tJig presence of neutron and proton cur- 
rents, the corresponding Fermi surfaces are displaced by 
a vector Q^. In the limit of small currents Qq ^ 
writing the "superfluid velocities" from as 



Vq = fiQ^/m 



(6) 



TO® = TO + &m®p + (5to®„ 



(14) 



where the various contributions to the effective masses 
are related to the Landau parameters by the simple for- 
mula M 



S-m-qq' = Z^^l ™W ® 

3 y TiqTO^ 



(15) 



The mass density matrix can then be equivalently writ- 
ten explicitly as 



TO -t- Smf 



-'nn — t'n 



(16) 
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Ppp — Pp' 



Pnp — Ppn 



8m® 



5mt 



Pp- 



(17) 



(18) 



It should be remarked that in the formulae provided by 
Sauls (see the terms proportional to 5ra®^ and (5m®, 
in the expressions for /5„„ and ppp are omitted and there- 
fore those formulae violate Galilean invariance. 

The quasiparticle energies e('^^{^} and the quasipar- 
ticle interaction {A;, A;'} can be deduced from a mi- 
croscopic approach. The solution of the many body 
problem, starting from the bare nucleon-nucleon inter- 
actions is very difficult. We shall here adopt a simpler 
approach based on self-consistent mean field models with 
phenomenological effective interactions (for a review, see 
for instance 10]), which have been very successful in de- 
scribing the nuclear properties of terrestrial nuclei. Such 
mean field models have also been widely applied in the 
context of neutron stars. 



III. LANDAU PARAMETERS IN THE ENERGY 
DENSITY FUNCTIONAL THEORY 

We shall calculate in this section the Landau param- 
eters for asymmetric nuclear matter in the framework 
of the Hohenberg-Kohn-Sham energy density functional 

theory [HIS- 

The energy density functional for spin-unpolarized ho- 
mogeneous nuclear matter is written as a sum of the 
isoscalar (T = 0) and isovector (T = 1) terms ^3 



T=0,1 



The isoscalar and isovector parts of some quantity for 
a nucleon system is given respectively by the sum and 
the difference between the neutron and proton contribu- 
tions. For example, the isoscalar and isovector densities 
are given by uq = Un + Up = Uh and ni = rin — Up 
respectively. 

The nucleon density n^, kinetic energy density Tq (in 
the units of Ti^ /2m), and nucleon current are express- 
ible in terms of the nucleon distribution function n^''^{k} 
by 



d^A: 

(27r)3 



(22) 



Energy density functionals of the form p9|l can be 
obtained in the Hartree-Fock approximation with effec- 
tive contact nucleon-nucleon interactions v{ri,r2\ of the 
Skyrme type, whose standard parametrisations (ignoring 
spin-orbit terms which are irrelevant in the present case) 
are 

v{ri,r2} = to(l + X(3Pa)5{ri -r2} 

+ iti(l + xiP^){k^^ 6{ri - ra} + 5{ri - r2]k) 

+ ^2(1 + X2Pa)k^ ■ S{ri - rajfc 
+ ^t,il+xsP.)S{r,-r2}n^{'^^±^}\ (23) 

where Pa — {1+ai •tT2)/2 is the spin exchange operator 
and k — — i(Vi — V2)/2. The density dependent term 
proportional to ^3 represents the effects of three body in- 
teractions. The coefficients C"{ri,f,}, and can then 
be expressed in terms of the parameters of the Skyrme in- 
teraction (see Appendix). It should be stressed however 
that the functional (|19|l is more general than the Skyrme 
functional. In particular, the coefficients C"{nb} can be 
any function of the baryon density n^,. 

The single particle energies are obtained from the func- 
tional derivative of the energy density 



{k} = 



S£ 



Snii){k} 



(24) 



where the zero subscript indicates that the functional 
derivative is evaluated in the static ground state (in 
which the currents vanish), characterized by the dis- 
tribution function 



,.0 >{k} = e{k'^'> - k} . 

Substituting the functional 119|) yields 



eq{k} 



2mf 



Ua 



(25) 



(26) 



in which the effective mass m® (using the same symbol 
as for the Landau effective mass defined by ^ since both 
definitions coincide) and the single-quasiparticle poten- 
tial Uq are given by 

68 



= — = — + (Q - Cl)n^ + 2Clnq (27) 



2to^ 5Tq 2m 



d^A: 

(2^ 



{k} 



d^A: 

(27r)3 



k^n^'i^k} 



(20) 



(21) 



U, = ^^ ^C^riq - 2nt{CS - C^) 

07lq 



(Q - CDrt + 2ClTq 



dnh 



dCf 
dnh 



{2nq 



(28) 
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where Tb = t„ + r^. jpp ^ 2k'i^\Cl + C[) + 2(Ci" + CJ) + ^4nb 

The quasiparticle interaction, calculated as a second drib 



functional derivative of the energy functional l^ni J is dC" d^CJ 2 d^Cf 

_ „ 



6ni<i){k}5n^i'){k'} 



{Un - Hp) + " "2" + " "2^ (»n - rip)^ , (31) 
d^b dn^ dnr 



(29) 



and contains only £ = and £ ^ I components. 

The non- vanishing Landau parameters are found to be 
expressible as 

/„"" = 24"^\C^ + CD + 2(Ci" + Co") + ^4nb 



+4^(n„ - Up) + ^^T%"b + ^^T%K - np)'^ , (30) 
cinh clrij^ arii^ 

I 



/7 = fP- ^ (4«)^ + fc(f)2)(CJ - CD + 2(Co" - CD + ^4nb + + ^("n " n^)^ (32) 



b dr.b 



/^ = 2(c•^ + c^)4")^ 



(P)2 



/r = 2(Ci5+cf)4' 



(33) 



(34) 



(35) 



These formulae agree with those of Bender et al. 
for the limiting case of symmetric nuclear matter (using 
standard notations = ff^ = fe+f'e and /"^ = ft—fe) 
and generalize the results of Blaizot & Haensel [Tjl for 
asymmetric nuclear matter to any energy density func- 
tional of the form p9l) . 

It should be remarked in particular that for any such 
functional (|19|l . the parameters /f" and /f^ are related 
by 



Jl 

/r 



2/3 



(36) 



The corresponding dimensionless £ = 1 Landau param- 
eters can be expressed in compact form as 



in which the coefficients aqqi are defined by 



n 



(37) 



(38) 



We conclude this section by remarking that in a general 
case of asymmetric nuclear matter (i.e., with n„ 7^ rip) 
the £ = 1 Landau parameters can be uniquely determined 
in terms solely of the effective masses as 



Up I 1 



m 



(40) 



m 

mi 



m 



(41) 



■pp _ 



n„ m 



m 



1 - 



(42) 



IV. ENTRAINMENT MATRIX AND 
EFFECTIVE MASSES 

Subtituting the expressions (|37|1 of the Landau param- 
eters obtained in the previous section, the entrainment 
matrix elements (|Hl can be seen to be expressible as 



Pqq' — Pq OA^qq' ^" '^qq'PqPq' i 

m„ 



(43) 



or more explicitly in terms of the mass densities p„ and 
Pp 

Pnn = Pn(l - a„pPp) (44) 



(39) 



Ppp — ^npPn) 



(45) 
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Pnp — ^npPnPp ■ 



(46) 



It is readily seen that the formulae (|44l4t)|) imply basic 
property of the entrainment matrix, namely p„„ + p„p = 
Pn and ppp + Pnp = Pp, which guarantees the Galilean 
invariance of the two fluid model. 

Two other kinds of effective masses, different from the 
Landau effective masses m® , have been introduced in the 
literature. Effective nucleon masses can be defined from 
the mass density matrix elements by setting 



Pqq 
Pq 



m 

™9 



(47) 



in such a way that in the proton momentum rest frame 
{V p = 0) we have 7r„ = ™j*Un and similarly in the neu- 
tron momentum rest frame, TTp = m^Vp. 

These effective masses have a very simple density de- 
pendence as shown on the formulae 



I 



m 1 - anpPp 



(48) 



These formulae show that, as pointed out by Prix et al. 
p^] . in the limit of very small proton fraction np/n\, <C 1, 
as relevant in the liquid core of neutron stars, we shall 
have m" ^ ~ m and m£ ~ mj*. We shall compute 
more accurate values of the effective masses in section 

ED 

In studies of neutron star cores, the non diago- 
nal entrainment matrix elements pnp has often been 
parametrised as 



Pnp 



-^Pn 



(55) 



in which the dimensionless parameter e was taken as a 
constant [H IHIll. Other authors [11 E i3| have 
suggested instead to set the dimensionless parameters de- 
fined by 



1 



a 

ml 



(56) 



as constants. However comparison with (I46II . H51|l and 
(|52|l shows that neither ansatz is satysfying, since these 
parameters are found to vary with the densities according 
to 



m 



(49) 



e — —OLnpPp ■ 



(57) 



The t|-effective masses differ from the Landau quasipar- 
ticle effective masses, and are related to the latter ones 
by 



OLnpPh 



(58) 



m 



m 



^qqPq ■ 



(50) 



due to the non-vanishing quasiparticle interactions. The 
extra term on the right hand side can be interpreted as 
resulting from the backflow induced by the motion of the 
quasiparticles. 

Alternatively one can introduce effective masses ml 
such that in the proton rest frame (meaning Vp = 0) we 
have TTn = m^Vn and similarly in the neutron rest frame 



These effective masses are given by 



n 



m 



p 

to; 



1 



^npPrt 



1 - OLnpPh 



1 



^npPp 



1 - anpPh 



(51) 



(52) 



where pb — Pn + Pp- The effective masses of the different 
kinds are related by 



Too- - m = (to" 



— m — (toJ — to) 



Up \ m 



(53) 



(54) 



^npPn 



1 



OLnpPh 



(59) 



The effective masses and entrainment parameters seem 
to diverge at some points of the p„ — pp plane. However 
as will be shown in the next section, once stability con- 
straints are imposed, these apparent singularities disap- 
pear. 



V. STABILITY OF THE STATIC GROUND 
STATE AND CONSTRAINTS ON THE 
ENTRAINMENT PARAMETERS 

Since the momentum of each nucleon is a linear com- 
bination of both the neutron and proton currents, this 
means that the corresponding dynamical contribution to 
the Lagrangian density of the system Adyn = i^^dyn is a 
bilinear symmetric form of the currents. It is therefore 
readily seen that this dynamical contribution is express- 
ible in terms of the Andreev-Bashkin entrainment matrix 
elements as 



1 2 

fdyn = l^iPnnV n 



2pnpVn ■ Vp + PppVt) . (60) 



As a result the total energy density £ of the fluid mixture 
can be written as the sum of the dynamical contribution. 
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£^dyn, and an internal static contribution Si^s which only 
depends on the densities: £ = fdyn + ^ins • 

The static ground state of the system is stable if the 
term Sdyn is strictly positive. This means that the en- 
trainment matrix must be positive definite (the minimum 
energy state thus being obtained by the vanishing of the 
superfluid velocities or equivalently of the currents, i.e. 
£dyn = 0), which mcaus that its eigenvalues must be 
strictly positive. This condition entails that the matrix 
elements (|44l46ll should obey 



Pnn ~i~ Ppp ^ ; 



Pnn ^ Pnn Ppp ■ 



(61) 



(62) 



These conditions lead to contraints on the £ = 1 Lan- 
dau parameters (using the other constraint that the Lan- 
dau effective masses m® and ra® have to be positive) 



(63) 



together with 



Whenever this condition is fulfilled, it can be seen from 
equations (|48I52|) that for any neutron and proton densi- 
ties, the effective masses ml and mj are therefore positive 
and smaller than the bare nucleon mass 

<ml,ml <m. (70) 

Combining the latter inequality with (|65|l shows in par- 
ticular that 



nq/n\y < ml/m < 1 . 



(71) 



Besides, since the proton fraction is very small inside 
neutron stars, it can be seen from H68|l and the definitions 
(|^ . 7 lEU and (|^ that in this case the neutron 
effective masses are always larger than the proton ones 



to" > ml , m^ > TOj . 



(72) 



Likewise, it can be shown that Sq > which, in associa- 
tion with (|67|l . yields 



< < 1 



(73) 



> 



(64) 



It is thus found that the entrainment parameters are 
well behaving functions of nucleon densities. 



The stability conditions can also be expressed in terms 
of effective masses i2L| 



ml n. 
> 



m 



q 



or equivalently 



TO, 
TO 



1 < !^ 



(65) 



(66) 



In terms of the dimensionless entrainment parameters Sq, 
these conditions can be expressed as 



Eg < 1 



iLq 

71b 



(67) 



It should be remarked that the previous inequalities are 
very general and have to be satisfied in any two fluid 
models. In the present case, these conditions also impose 
a constraint on the energy functional (|19|) from which the 
entrainment matrix is derived. Since the conditions H61|) 
and (|52|l must be satisfied for any neutron and proton 
densities, this leads to the following requirement 



(68) 



In the particular case of Skyrme functionals, this last 
condition reads 



VI. APPLICATION TO NEUTRON STAR 
MATTER IN BETA EQUILIBRIUM 

In the previous section we have obtained general for- 
mulae for the entrainment parameters and the associated 
effective masses for nuclear matter with arbitrary asym- 
metry. In the present section we shall apply these for- 
mulae to construct a model of neutron star core. We 
shall use the SLy4 Skyrme force which has been specifi- 
cally devised for astrophysical purposes ,22c 2i 2^ 2.3 ■ 
In the framework of a compressible liquid drop model 
based on the SLy4 Skyrme energy functional, Douchin 
& Haensel 26J found that the bottom edge of the crust 
corresponds to the baryon density riodge — 0.076 fm~'^. 
In the following we shall consider the density domain 
"■edge < 'Jb < 3ns, where — 0.16 fm^"^ is the nuclear 
saturation density. 

We assume that the liquid core is composed of an ho- 
mogeneous plasma of neutrons, protons and electrons 
(and muons for baryon densities beyond some critical 
threshold) in beta equilibrium 



(74) 



This means that the chemical potentials of the various 
species are related by (assuming that neutrinos escaped 
from the star) 



ti{2 + Xi) + t2{2 + X2) > 



(69) 



f^n — ^~ Me ; Me — Mm 



(75) 
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In the Hartree-Fock approximation, chemical potentials 
of nucleons are equal to the corresponding Fermi energies 
{q — n,p) including rest mass energy 

flq = rUqC^ + ^ + Uq . (76) 

ziriq 

Considering the leptons as ideal relativistic Fermi gases, 
the lepton chemical potentials are given by {I — e, ^) 

^ y^m,c2 + ;i2^2(37,2„^)2/3. (77) 

Charge neutrality requires that 

Up ^ He + . (78) 
I 



In equations (|76|l and (|77|l . we have neglected the de- 
viations in the chemical potentials due to the existence 
of non vanishing currents since the relative velocities are 
typically very small compared to the velocities of the var- 
ious constituents. For completeness, let us mention that 
the internal energy density of the nucleons can be decom- 
posed in the form 



^int — ^^0 + ^cnt (79) 

in which £q is the functional (|19|l evaluated in the static 
ground state with the distribution function H25|l including 
the rest mass energies. 



, m 

I 



and Sent is the entrainment contribution which is express- 
ible as 

font = ~^Pn£n{Svf (81) 

where Sv is the velocity difference between neutrons and 
protons. The entrainment term is negligibly small com- 
pared to the static term E^nt ^ £o even for the fastest 
pulsars and can therefore be neglected. 

The muons are present in matter when the electron 
chemical potential fie exceeds the muon mass m^c^ ~ 
105 MeV. This occurs at a baryon density rib — 0.12 
fm^''. In equilibrium, the composition of the liquid core 
is therefore completely determined by the baryon density 
Ub and is shown on figure ^ 

The dimensionless entrainment parameters as defined 
by H55|l and H56|) . which have been widely used in neu- 
tron star simulations, are represented on figures [5] and |31 
respectively. The and [I- effective masses are shown on 
figures 0] and [51 respectively. Due to the increase of the 
proton fraction with the baryon density (see figure^ , the 
differences between the two kinds of effective masses ml 
and TO J, which are negligible at the crust-core boundary, 
become significant in deeper layers. We have also plotted 
on figure IHl the Landau effective masses for comparison. 
As can be seen on those figures, the various definitions 
of "effective mass" do not coincide. This concept should 
therefore be carefully employed and the definition that 
has been adopted should always be clearly specified. 

We have finally shown on figure [3 the dimensionless 
determinant of the entrainment matrix 

2 

J. _ PnnPpp — Pnp ^-g^-j 
PnPp 



which appears in the perturbed hydrodynamical equa- 
tions and which is therefore important for the study of 
oscillation modes. In the present case, this quantity is 
simply given by 

T = 1 - UnpPb , (83) 

and is therefore quite remarkably independent of the nu- 
clear asymmetry. 

For comparison, we have computed the coefficient 
cinp from which all the entrainment parameters can be 
obtained, for the 27 Skyrme forces recommended by 
Rikovska Stone et al. for neutron star studies. 

The coefficient 5„p ranges from for the parametrisa- 
tions SkT4 and SkT5 ^ down to -10.4168 (in units 
m^ ^.fm^^) for the parametrisation SV J29|. Let us also 
mention that the Skyrme SLy forces iH El HI El, 
which have been widely employed in neutron star stud- 
ies, yield coefficients around (5„p ~ —1.5 except for the 
SLy230a force for which 5„p = -0.007360. 

Let us notice that these three forces SkT4, SkT5 and 
SLy230a are the only parametrisations for which the 
isovector effective mass, relevant for the T = 1 isovec- 
tor electric dipole sum rule, was set equal to the bare 
nucleon mass (see Table VII of reference H^)- However, 
the isovector electric dipole "giant resonance" in nuclei 
consists essentially of relative motion of protons against 
neutrons, and the sum rule constraint is therefore crucial 
for the entrainment effect in the infinite neutron-proton 
mixture. The case a„p = implies that the effective 
masses to^ and ml are equal to the bare one and there- 
fore there is no entrainment. Reciprocally large negative 
values of associated with strong entrainment ef- 

fects. 
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FIG. 1: (Color online) Equilibrium fractions of protons 
(dashed line), electrons (thick line) and muons (thin line) in 
neutron star liquid core as a function of the baryon density 
rib = rip + n„ from the bottom edge of the crust ricdgc — 0.076 
fm""^ down to Sris where rig = 0.16 fm"'' is the nuclear satura- 
tion density. The results have been obtained with the Skyrme 
SLy4 energy density functional. 
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FIG. 3: (Color online) Dimensionless entrainment parame- 
ters En (solid line) and Sp (dashed line) for npe/i matter in /3 
equilibrium (Skyrme SLy4 energy functional). 
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FIG. 4: (Color online) Effective masses m"/m (solid line) and 
mZ/m (dashed line) for npen matter in /3 equilibrium (Skyrme 
SLy4 energy functional). 



FIG. 2: (Color online) Dimensionless entrainment parame- 
ter e for npefi matter in /3 equilibrium (Skyrme SLy4 energy 
functional). 



VII. CONCLUSION 



In the present paper we considered only spin- 
unpolarized nuclear matter. We therefore did not dis- 
cuss spin and spin-isospin instabilities that plague many 
Skyrme forces at supranuclear densities [30, • How- 
ever, let us notice that for the SLy4 used in our calcu- 
lations, the ferromagnetic instability appears above the 
baryon density 0.5 fm^'^, which is beyond the upper limit 
in our figures. 



Analytical expressions for the entrainment matrix and 
related effective masses of a neutron-proton superfluid 
mixture at zero temperature have been obtained within 
the non relativistic energy density functional theory. In 
contrast to recent investigations within relativistic mean 
field models 4], the entrainment parameters have been 
found to be expressible by very simple formulae which 
could be easily implemented in dynamical simulations 
of neutron star cores. We have also clarified the link 
between various definitions of effective masses that have 
been introduced in the literature. 

We have applied these formulae for Skyrme forces in 
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FIG. 5: (Color online) Effective masses m^/m (solid line) 
and m^/m (dashed line) for npe/i matter in (3 equilibrium 
(Skyrme SLy4 energy functional). 



APPENDIX A: SKYRME ENERGY DENSITY 
FUNCTIONAL COEFFICIENTS 



The energy functional deduced from the Skyrme effec- 
tive interaction in the Hartree-Fock approximation has 
a similar form as equation l|19|l . The coefficients in the 
energy functional (|19() can thus be expressed in terms of 
the parameters of the Skyrme interaction as follows. As 
a result of local phase invariance of the Skyrme forces 
[s^. the coefficients C;j, and are related by 



Crp Crp . (Al) 



In terms of the parameters of the Skyrme interaction the 
coefficients of the energy functional are given by 




FIG. 6: (Color online) Landau effective masses m®/m (solid 
line) and mf /m (dashed line) for npe/x matter in (3 equilib- 
rium (Skyrme SLy4 energy functional). 



order to evaluate the entrainment matrix in the stan- ^ ''2 ~^ '''■^'^ ~^ 8^^'' 2 (^^) 

dard model of the liquid core of neutron stars, composed 
of a mixture of neutrons, protons, electrons and possi- 
bly muons in beta equilibrium. Comparing the results 
with different Skyrme forces, we have found that the en- 
trainment parameters are quite sensitive to the adopted 
parametrisation. The observations of the entrainment 

effects in neutron stars could therefore provide new con- ACKNOWLEDGMENTS 
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